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ABSTRACT 

We measure the morphology-density relation (MDR) and morphology-radius relation (MRR) for 
galaxies in seven z ^ I clusters that have been observed with the Advanced Camera for Surveys on 
board the Hubble Space Telescope. Simulations and independent comparisons of our visually derived 
morphologies indicate that ACS allows one to distinguish between E, SO, and spiral morphologies 
down to Z850 = 24, corresponding to L/L* = 0.21 and 0.30 at z = 0.83 and z = 1.24, respectively. We 
adopt density and radius estimation methods that match those used at lower redshift in order to study 
the evolution of the MDR and MRR. We detect a change in the MDR between 0.8 < z < 1.2 and 
that observed at z ~ 0, consistent with recent work - specifically, the growth in the bulge-dominated 
galaxy fraction, /e+so; with increasing density proceeds less rapidly at z 1 than it does at z ~ 0. 
At z 1 and S > 500 galaxies Mpc"^, we find </e-hso>= 0.72 ±0.10. At z - 0, an E-fSO population 
fraction of this magnitude occurs at densities about 5 times smaller. The evolution in the MDR is 
confined to densities E > 40 galaxies Mpc~^ and appears to be primarily due to a deficit of SO galaxies 
and an excess of Spiral-|-Irr galaxies relative to the local galaxy population. The /e - density relation 
exhibits no significant evolution between z = 1 and z = 0. We find mild evidence to suggest that 
the MDR is dependent on the bolometric X-ray luminosity of the intracluster medium. Implications 
for the evolution of the disk galaxy population in dense regions are discussed in the context of these 
observations. 

Subject headings: galaxies: clusters: general — galaxies: formation — galaxies: evolution — galaxies: 
structure 



1. INTRODUCTION 

The study of the origin and evolution of the morpho- 
logical distribution of galaxies in different environments 
can reveal important information about internal galactic 
stellar and gas dynamics, the state of star formation ac- 
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tivity as a function of time, and constrain the relative sig- 
nificance of the effects of environmental processes versus 
conditions at the epoch of their formation in establish- 
ing galactic structure. In standard hierarchical clustering 
models, galaxies in high density regions of the Universe, 
such as in the central regions of galaxy clusters, will col- 
lapse earlier and may evolve more rapidly than galax- 
ies in low density regions (Kauffmann 1995; Benson et 
al. 2001, Heavens et al. 2004). In addition, galaxies in 
dense environments are subject to a variety of external 
stresses, which are, in general, not conducive to the main- 
tenance of spiral structure. These processes include ram- 
pressure stripping of gas (Gunn & Gott 1972; Farouki 
& Shapiro 1980; Kent 1981; Fujita & Nagashima 1999; 
Abadi, Moore Sz Bower 1999; Quilis, Moore & Bower 
2000), galaxy harassment via high speed impulsive en- 
counters (Moore et al. 1996, 1999; Fujita 1998), cluster 
tidal forces (Byrd & Valtonen 1990; Valluri 1993; Fujita 
1998) which distort galaxies as they come close to the 
center, interaction/merging (Icke 1985; Lavery & Henry 
1988, Mamon 1992; Makino & Hut 1997; Bekki 1998), 
and removal and consumption of the gas due to the clus- 
ter environment (Larson, Tinsley & Caldwell 1980; Bekki 
et al. 2002). 

Two of key relationships that must be understood in 
the context of the above processes are the relative pop- 
ulation fraction of the different morphological classes 
as functions of the local galaxy density and their loca- 
tion within the local gravitational potential well. The 
morphology - density relation (hereafter MDR) and 
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the morphology ~ radius relation (hereafter MRR) have 
been well studied at low-z (Dressier 1980 - hereafter 
D80; Postman & Geller 1984 - hereafter PG84; Whft- 
more & Gilmore 1991; Goto et al. 2003a) and quantify 
many long-standing observations showing a preference 
for spheroidal systems to reside in dense regions (or per- 
haps better stated as a significant lack of spiral galaxies 
in dense regions). A full understanding of how such a 
cosmic arrangement came to be requires measuring the 
evolution of the MDR and MRR. Such an evolutionary 
study is only possible using the high angular resolution 
provided by the Hubble Space Telescope (HST). Several 
pioneering works have now shed light on this evolution 
(Dressier et al. 1997 - hereafter D97; Fasano et al. 2000; 
Treu et al. 2003; Smith et al. 2004 - hereafter Sm04). 
D97 and Fasano et al. (2000) find a significant decline in 
the fraction of lenticular galaxies (/so) when one looks 
back from the present epoch to an epoch 4—5 Gyr ago. 
The results presented by Treu et al. 2003 and Sm04 are 
perhaps the most enlightening - they find a smaller in- 
crease in the bulge-dominated galaxy (E-I-SO) fraction 
(/e+so) with increasing density at z > 0.4 than is seen 
at z < 0.1 but also find comparable /e+so values for 
low-density regions (S < 10 galaxies Mpc"^) at z > 0.4 
and the current epoch. Sm04 propose a simple model to 
explain these observations in which high density regions 
at z ~ 1 would largely be comprised of elliptical galaxies 
with only a trace of lenticulars {e.g., </so< O-l)- They 
consider various processes to transform spiral galaxies 
into lenticular galaxies in order to increase /so with time 
to match the observed morphological population frac- 
tions at 2 ~ 0.5. However, the Sm04 z ~ 1 /so mea- 
surement was inferred from /e+so rather than measured 
directly as Sm04 chose (perhaps wisely) not to attempt to 
distinguish between SO and E galaxies from the WFPC2 
data used in their study. 

The deployment of the Advanced Camera for Surveys 
(ACS; Ford et al. 2003) on the HST has provided us 
with an opportunity to greatly expand our understand- 
ing of the physics behind the morphological evolution of 
galaxies in a wide variety of environments. The higher 
sensitivity and better angular sampling of the Wide Field 
Camera (WFC) on ACS relative to WFPC2 enables the 
acquisition of high S/N morphological information for 
sub-i* galaxies over projected areas of up to 10 Mpc^ in 
z '--^ 1 clusters, in a modest allocation of telescope time. 
This is a significant advantage over prior capabilities and 
enables us to sample > 3 decades in local galaxy density 
using the same homogeneous data samples. 

As part of an extensive program to study the formation 
and evolution of clusters and their galaxy populations, 
the ACS Investigation Definition Team (IDT) has imple- 
mented a 128 orbit program to observe 7 distant clusters 
in the redshift range 0.83 < z < 1.27. In this paper, 
wc present new constraints on the form and evolution of 
both the MDR and the MRR in these clusters and their 
surroundings based on morphologies determined from the 
ACS/ WFC imagery coupled with extensive spectroscopic 
data and X-ray observations. This paper is organized as 
follows: section 12 contains a brief summary of the space 
and ground-based observations used in this study, sec- 
tioninipresents a detailed discussion of our morphological 
classification procedure and an assessment of the reliabil- 
ity of these classifications, section0]presents the methods 



used to estimate the local projected density, and section[Sl 
presents our derived MDR and MRR. An assessment of 
the implications of these results is given in section|Hland a 
summary of the essential results is provided in section |7| 
Two appendices discuss details associated with the com- 
putation population fractions that are suitably corrected 
for contamination and incompleteness, the robustness of 
our density estimators, and the validity of using com- 
posite samples to enhance the signal-to-noise ratio in the 
derived MDR and MRR. We adopt Ho = 70 km s'^ 
Mpc^^, n„i — 0.3, and Ha — 0.7 for the computation of 
all intrinsic quantities unless specifically indicated other- 
wise. 

2. OBSERVATIONS 

The clusters included in this study, along with a sum- 
mary of the ACS observations, are listed in Tabled The 
average cluster redshift, based on all available spectro- 
scopically confirmed cluster members, is given in column 
2 of this table. The number of redshifts acquired for 
galaxies in the region of each cluster (both members and 
non-members and including galaxies outside the ACS 
mosaic boundaries) is listed in column 3. Column 4 lists 
the number of spectroscopically confirmed cluster mem- 
bers that also lie within the ACS mosaic boundaries. The 
details of the HST ACS observations are given in columns 
5-7. The sample selection process was limited by the 
small number of spectroscopically confirmed clusters at 
z > 0.8. However, we were able to include clusters with 
a range of X-ray luminosities, from Lx.Bo1< 10^* erg s~^ 
to ix,Boi= 1-9 X lO''^ erg s~^. Table|21is a compilation 
of the derived X-ray properties and velocity dispersions 
of these clusters. Two of the seven clusters (the two at 
R.A. ~ 16 hr) are optically selected systems, the rest 
are X-ray selected although RXJ0849-I-4452 is a binary 
cluster system in which the less massive component (CL 
J0848-H4453) was IR-selected (Stanford et al. 1997). 

2.1. ACS Observations 

We used the WFC on the ACS to image each clus- 
ter. For MS1054-0321, RXJ0152-1357, RDCS1252-2927, 
and RXJ0849-)-4452 multiple pointings were used to con- 
struct larger mosaics covering ~ 35 square arcminutes. 
For the first three of these clusters, the pointings form a 
2x2 pattern with all 4 pointings overlapping the central 
1 arcminute region of each cluster - hence the exposure 
time in the central regions of these systems is 4 times as 
long as the values given in TableHl For RXJ0849-t-4452, 
we used a 3 x 1 pattern in order to obtain images of both 
components of this binary cluster system. All the re- 
maining targets were observed using a single WFC point- 
ing centered on the cluster. 

The filters are chosen to approximately straddle the 
rest-frame 4000A break in order to facilitate the identifi- 
cation of bulge-dominated galaxies in the red sequence of 
the clusters. This sequence, which is populated mostly 
by elliptical and lenticular galaxies with a strong 4000A 
break, is well separated from the CM relation for late- 
type cluster galaxies as well as that for most field galax- 
ies. In all cases, we have at least one filter that samples 
part of the rest-frame i?— band. We use the ACS im- 
ages taken in those filters to perform our morphological 
classifications so that we can readily compare our re- 
sults with morphological information obtained at lower 
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redshifts (e.g., Fabricant, Franx, & van Dokkum 2000). 
The filters used are explicitly listed in Table H Here- 
after, we will use Vqoq to denote the F606W bandpass, 
''625 to denote F625W, 1775 to denote F775W, Isu to 
denote F814W, and zsso to denote F850LP. 

2.2. Object Photometry and Classification 

Object detection and analysis is performed using the 
ACS IDT Pipeline {a.k.a. APSIS; Blakeslee et al. 2003a). 
APSIS photometry is on the AB system and is corrected 
for Galactic extinction using the Schlegel et al. (1998) 
lOOfi map. Object detection and star-galaxy discrimina- 
tion is done using the dual image mode in SExtractor 
(Bertin & Arnouts 1996). The detection image is an 
inverse variance weighted combination of the ACS ex- 
posures from all available passbands. The inverse vari- 
ance weighting preserves information about the struc- 
tural characteristics of the galaxies (see Benitez et al. 
2004 for details). In this paper, we count as galaxies all 
objects with SExtractor stellarity parameter CLASS_STAR 
< 0.50. The automated image structure analysis of de- 
tected objects in our ACS data includes the determina- 
tion of the luminosity-weighted moments, the ellipticity, 
and the 180° rotational asymmetry and image concen- 
tration parameters {e.g., Abraham et al. 1994; Conselice 
et al. 2000). All magnitudes cited in this study are based 
on the SExtractor MAG_AUTO magnitude as it provides a 
reasonable estimate of the total flux. 

2.3. Spectroscopic Observations and Photometric 
Redshifts 

Spectroscopic redshifts have been obtained for the clus- 
ters in our survey, by us and others, using multiob- 
ject spectrographs on the Keck, VLT, or the Magel- 
lan observatories. The total number of redshifts and 
the number of confirmed cluster members within each 
ACS mosaic are listed in Tabled The publications con- 
taining some or all of the redshift data include Tran 
et al. (2005) for MS1054-0321, Demarco et al. (2004) 
for RXJ0152-1357, Postman et al. (1998a, 2001) and 
Gal & Lubin (2004) for the CL1604-I-43 system, Stan- 
ford et al. (2002) for RDCS0910-H5422, Demarco et al. 
(2005) for RDCS1252-2927, and Rosati et al. (1999) for 
RXJ0849-I-4452. The target selection criteria for the red- 
shift surveys of MS1054-0321 and the CL1604-H43 clus- 
ters were based on a single red flux limit. The target 
selection for RXJ0152-1357 include a color-selection cri- 
terion (see Demarco et al. 2004). The spectra are of 
moderate resolution (R ~ 500 — 1200) and most have suf- 
ficient S/N to measure the prominent spectral features 
(e.g., [Oil] line widths). 

Our photometric redshifts are derived using the 
Bayesian method {a.k.a. BPZ) described in Benitez 
(2000) and are based on a minimum of 3 passbands, in- 
cluding all available ACS photometry. We have reliable 
photo-z's for RXJ0152-1357 {z = 0.837), MS1054-0321 
{z = 0.831), and RDCS1252-2927 {z = 1.235). Photo-z's 
for RXJ0152-1357 are based on the ACS r625, 1775, and 
Z850 photometry. Photo-z's for MS1054-0321 are based 
on the ACS Vgoei *775, and zgso photometry. Photo-z's 
for RDCS1252-2927 are based on ACS ^775, zgso pho- 
tometry and ground-based BVJK photometry (Toft et 
al. 2004). For the two z = 0.83 clusters, the rms scatter 



in {zspec-Zph)/{l + Zspec), cTph, is 0.05. For RDCS1252- 
2927, aph is 0.10. Only objects that have a BPZ con- 
fidence level of 0.90 or greater are selected for analysis. 
Figure n shows the distribution of galaxies with photo- 
metric redshifts within 2(Tp/i of the mean cluster redshift 
and galaxies with 2aph < \zcl - Zph\/{1 + zcl) < GcFph- 
The cluster overdensity is clearly seen only when we se- 
lect galaxies close to the actual mean cluster redshift, in- 
dicating our photo-z's are useful in significantly suppress- 
ing fore/background contamination and isolating most of 
the actual cluster members. 

3. MORPHOLOGICAL CLASSIFICATION 

We visually classified the morphologies of all galaxies 
in each field with 1775 < 23.5 (for the z < 1 clusters) 
or Z850 < 24 (for the z > 1 clusters) regardless of their 
position or color. For reference, the characteristic mag- 
nitude, m*, for cluster galaxies is 1775 = 22.3 at z — 0.83 
(Goto et al. 2004) and zggo = 22.7 at z = 1.24 (Blakeslee 
et al. 2003b). For all our cluster observations, we have 
at least one filter that samples part of the rest-frame 
i3— band (see Table ^ so that morphological classifica- 
tions can be readily compared with those at lower red- 
shifts. We classify galaxies using the common Hubble 
sequence: E, E/SO, SO, SO/Sa, Sa, Sb, Sb/Sc, Sc/Sd, 
Irr. However, for the purposes of the present analyses, 
we bin these finer classifications into just 3 broad mor- 
phological categories: E (elhptical; — 5 < T < —3), SO 
(lenticular; — 2 < T < 0), and Sp-|-Irr (Spiral -t- Irregular; 
1 < T < 10). The FWHM of the point spread function 
in our co-added ACS images is ~ 0.09 arcsec (1.8 WFC 
pixels), corresponding to a projected proper distance of 
684 pc at z = 0.831 and 752 pc at z = 1.27. We are thus 
able to resolve sub-kpc structure in all cluster members. 
At 1775 = 23.5, the typical galaxy subtends an isophotal 
area of -400 WFC pixels or 125 x (FWHM)2, making 
visual classification (or for that matter any classification 
method) quite feasible. At fainter magnitudes, however, 
classification rapidly becomes difficult and systematic er- 
rors increase both because galaxies are becoming smaller 
(e.g., Roche et al. 1998; Bouwens et al. 1998; Ferguson 
et al. 2004; Trujillo et al. 2004) and because there is in- 
sufficient area over which the integrated S/N is sufficient 
for accurate classification. Examples of the ACS image 
quality and our corresponding classifications are shown 
in Figures 121 and 131 

The morphological classification was performed on the 
full sample of 4750 galaxies (in 7 clusters) by one of us 
(MP). Three other team members (NC, MF, BH) clas- 
sified a subset of ^ 400 of these galaxies to provide an 
estimate of the uncertainty in the classifications. All clas- 
sifiers used a common reference set of morphologies from 
a low redshift i?— band galaxy sample as a guide. Ex- 
act or majority agreement between all 4 classifiers in the 
overlap sample was typically achieved for 75% of the ob- 
jects brighter than 1775 = 23.5. Furthermore, there was 
no significant systematic offset between the mean classi- 
fication for the 3 independent classifiers (as determined 
using the voting scheme from Fabricant et al. 2000) and 
the classification by MP giving confidence that the full 
sample was classified in a consistent manner. 

The overall population fractions between the 4 inde- 
pendent classifiers exhibit only a relatively small vari- 
ance. Figure 21 shows the E-fSO fraction for each classi- 
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Fig. 1. — The projected distribution of galaxies with photometric redshifts that he within it2a-ph of the mean cluster redshift are shown 
in the right-hand panels for RXJ0152-1357, MS1054-0321, and RDCS1252-2927. The distribution of galaxies with photometric redshifts 
in the range 2crp^ < I^CL — ^ph\/{^ + -^Cl) < ^cTph are shown in the left-hand panels. Different symbols denote different morphological 
classifications: black dots are ellipticals, grey dots are SO's, and stars are Sp+frr. The dashed lines denote the boundaries of the ACS 
mosaics. The RDCS1252-2927 photo-z's are available over less area than the full ACS mosaic because they rely on NIR photometry that 
covers a smaller region. The overdensities associated with the clusters are easily seen in the right-hand panels and arc dominated by E and 
SO galaxies. 



fier as a function of Zgso niagnitude for a color-selected 
((*775 - Zs5o) > 0.5) subset of galaxies in the RDCS1252- 
2927 field. The average rms scatter in the E-I-SO pop- 
ulation fraction between classifiers is 0.06. The average 
rms scatter in the SO population fraction between clas- 
sifiers is about 2 times higher, 0.11. The rms in the E 
population fraction is the same as that for the SO's, 0.11. 
In other words, the population fractions are fairly robust 
and the variance in these fractions is significantly less 
than the ^ 20 — 25% disagreement level between classi- 
fiers on the morphological classification of any individual 
galaxy. Not surprisingly, the combined E-I-SO fraction 
is more robust than either the E or the SO fractions ~ 
a quantitative demonstration that detecting spiral struc- 
ture is a more robust skill than detecting disks. Most im- 



portantly, there are no significant systematic differences 
between the classifiers. The scatter in the SO popula- 
tion fraction is comparable with ^/N uncertainties in any 
given value and while counting statistics do not suggest 
the minimum scatter one might expect between different 
classifiers, this level of scatter indicates that our classi- 
fication errors are small enough for the task at hand - 
providing a uniform set of classifications. 

If our visual morphological classifications are robust, 
there should be noticeable differences in the distributions 
of the objectively derived "form" parameters (ellipticity, 
asymmetry, concentration) for the E, SO, and Sp-|-Irr cat- 
egories. Figure [S] shows the histograms of these 3 form 
parameters for each of the 3 morphological bins. Clear 
differences between the distributions exist. For example, 
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Full Resolution Figure Available at 
http://www.stsci.edu/~postman/mdr_figures 



Fig. 2. — Color postage stamp cutouts of the brightest 49 spectroscopically confirmed members of CL1604+4321 (z = 0.92) and 
CL1604+4304 (z = 0.90) and their corrcsijoiiding visually derived morphological classifications. The "Pos Disk" classification stands 
for "possible disk" galaxy. It is given to objects that appear to have a disk structure but the precise nature of that disk could not 
established. Galaxies with the "Pos Disk" classification are counted as Spirals in our derivation of the MDR. The first 29 cells show 
galaxies from CL1604+4321 and the rest contain galaxies from CL1604+4304. The galaxies from each cluster are shown in increasing 
apparent magnitude order and span the range 20.91 < 1775 < 23.50. Each cutout subtends a 6.4 x 6.4 arcsecond area. 



the ellipticity distribution of visually classified elliptical 
galaxies differs from that for visually classified SO galax- 
ies at greater than the 99.999% confidence level. We will 
use this fact later on as a key part of our analysis. Ellip- 
ticals as a class have, as expected, a lower median ellip- 
ticity and asymmetry and a higher median concentration 
than the other two morphological classes. Ellipticals also 
exhibit less scatter about the mean values of these form 
parameters. The Sp-|-Irr class exhibits a higher mean el- 
lipticity and asymmetry than either the E or SO class. 
The SO galaxies have form characteristics that arc, on 
average, intermediate between the E and Sp+Irr distri- 
butions. 

Two key concerns when performing galaxy morpho- 
logical classification, visually or via machine algorithms, 
are the effect of surface brightness dimming and the 



shorter rest-frame wavelength being imaged with increas- 
ing redshift. The latter effect, sometimes referred to as 
the "morphological k-correction," has been studied fairly 
well (Bunker et al. 2000, Abraham & van den Bergh 2001, 
Windhorst et al. 2002, Papovich et al. 2003). In general, 
over the redshift range being studied here, the morpho- 
logical k-correction has been shown to be important only 
in a small fraction 20%) of the galaxy population and 
in those cases, it is often more of an issue of how one char- 
acterizes any existing spiral structure and not usually a 
case of missing spiral structure altogether (sec references 
above for details). Furthermore, the bluest wavelengths 
we use for our morphological classifications correspond to 
the blue end of the rest-frame i?— band, which is where 
many lower redshift studies have been done. Our classi- 
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Full Resolution Figure Available at 

http:/ /www. stsci.edu/~postman/mdr_figures 



Fig. 3. — Color postage stamp cutouts of the brightest 49 spectroscopically confirmed members of MS1054-0321 (2 = 0.831). The 
galaxies are displayed in increasing apparent magnitude order and span the range 20.14 < 1775 < 22.15. 



fications are never performed in the rest-frame t/— band. 

The effects of SB dimming are potentially of greater 
concern as the ability to distinguish between adjacent 
categories {e.g., E vs SO, SO vs Sa) may be compromised 
at higher redshift. To test how sensitive our visual clas- 
sification scheme is to SB-dimming, we performed two 
simulations. In the first test, we "redshift" our ACS im- 
age of the z = 0.33 cluster MSI358-I-6245 to z = 0.83 and 
perform visual classifications on both the original and 
redshifted versions. In the second test, we redshifted our 
ACS image of MS1054-0321 to z = 1.24 and compared 
the classifications derived for the original and redshifted 
versions. The "redshifting" involved dimming the ob- 
jects appropriately, resampling the images to account for 
the smaller angular scale, and adjusting the noise lev- 
els to correspond to those appropriate for our exposure 
times used in the more distant cluster observations. Re- 
classification of the redshifted galaxy images was done 
in a random order and at least 3 months after the initial 



classifications to minimize "memory" of the initial classi- 
fications by the classifier. Some examples of the original 
and redshifted MS1054-0321 galaxies are shown in Fig- 
ure El The results of the comparisons between the mor- 
phological classifications of the original and redshifted 
objects are shown in Figure [7| The population fractions 
obtained using the redshifted images are completely con- 
sistent with those in the original images - the differences 
are comparable to or less than the ^/N uncertainties. 
Thus, our increased exposure times for the more distant 
clusters coupled with the high angular sampling and sen- 
sitivity of the ACS/WFC successfully mitigate the effects 
of SB-dimming and allow us to distinguish between our 
3 primary morphological categories (E, SO, and Sp-|-Irr) 
uniformly across the redshift range under study. 

3.1. Classification of Mergers & External Morphology 
Comparisons 
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Fig. 4. — Upper panel: The E+SO population fraction as a function of 2350 magnitude for each of the 4 independent visual classifiers. 
The sample used in this comparison is a color-selected sample that favors inclusion of bulge-dominated galaxies, hence the relatively high 
overall early-type population fraction. The values for each classifier are offset slightly from one another for clarity. The error bars shown 
represent the \/~N uncertainties. The dashed line shows the mean E-I-SO population fraction averaged over all classifiers. Lower panel: 
The RMS scatter in the E-I-SO (stars) and SO (diamonds) population fractions as a function of zsso magnitude between the 4 classifiers. 
The dashed lines represent the average RMS scatter. 



As we wish to measure the MDR and MRR in forms 
that can be compared to previous work, we attempt to 
provide a standard Hubble type classification (E, SO, or 
Sp+Ir) for all galaxies above our flux limits. We do 
make a separate note on whether the object appears 
to be undergoing a merger or tidal disruption and the 
analyses of the distribution and frequency of such sys- 
tems will be presented in a separate paper (Bartko et al. 
2005). However, for the present work, we do not clas- 
sify galaxies as merger /peculiar systems (as done by van 
Dokkum et al. 2000; hereafter referred to as vD2000) if 
one of the above standard morphological categories can 
indeed be applied to the individual objects involved in 
the merger. Nonetheless, the WFPC2 study of MS1054- 
0321 by vD2000 provides an additional check on the ro- 



bustness of our morphological classifications. There are 
a total of 79 galaxies in common that have morphological 
classifications by us and by vD2000. Of these, 16 are clas- 
sified as M/P (merger/peculiar) by vD2000. We classify 
13 of them as bulge-dominated systems (E, SO, or SO/Sa) 
and 3 as later- type spiral galaxies. The ACS cutouts of 
these 16 objects are displayed in Figure |H1 As can be 
seen from this figure, the morphology of most of the sys- 
tems classified as mergers/peculiar by vD2000 can also 
reasonably be placed into one of the E/SO/Sp-|-Irr bins. 
While this confirms the vD2000 conclusion that MS1054- 
0321 hosts a significant fraction of early-type mergers, it 
does yield one difference (albeit perhaps a semantic one) 
in the conclusions reached regarding the overall early- 
type populations in MS1054-0321. By counting merg- 
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Fig. 5. — The image concentration, rotational asymmetry, and ellipticity distributions for galaxies visually classified as elliptical, SO, 
and Sp+Irr. The top row shows the distributions for Sp+Irr. The central row shows the SO distributions. The bottom row shows the 
distributions for ellipticals. 



ers/peculiar systems as a separate category, vD2000 con- 
cluded that early-type systems comprise a lower popula- 
tion fraction (44%) than in comparably rich clusters at 
lower redshift. We conclude that the early-type fraction, 
/e+soj ill MS1054-0321 is higher about 73%, when one 
attempts to classify cluster members (including merger 
components) as E, SO, or Sp-|-Irr. Of course, the factor 
is a function of local density and the above value is aver- 
aged over densities in the range 15 < S < 1000 galaxies 
Mpc^^. Given that the merger/peculiar category has not 
routinely been used in the classification of low- 2; clusters, 
we feel it is an important exercise to provide a set of mor- 
phological classifications that are as similar as possible to 
the \ow-z studies. We also agree, however, that quanti- 
fying the frequency of mergers at higher redshift reveals 
fundamental information about the evolution of cluster 
galax;ies. The assessment of the early-type population 
component of MS1054-0321 by us and vD2000 are ulti- 



mately consistent, however, if one accounts for the ob- 
servation that the majority of the close pairs in MS1054- 
0321 include at least one bulge-dominated member. 

We compare our Hubble classifications with those 
galaxies that vD2000 did classify as E, SO, or Spiral/Irr 
(i.e., excluding the M/P objects). We find exact agree- 
ment with their Hubble classification (when binned into 
these 3 categories) 71% of the time. We swapped E or 
SO classifications 11% of the time (i.e., they called it E 
and we called it SO or vice versa) and wo swapped SO and 
Spiral classifications 10% of the time. The remaining 8% 
were cases where either we or they could not make a 
definitive classification. This translates to a ±0.1 scatter 
between our respective /e or /so values, which is con- 
sistent with the scatter estimated from our comparisons 
between our ACS team classifiers. 

As a further external validation of our morphologi- 
cal classifications, MP classified all galaxies from our 
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Full Resolution Figure Available at 

http:/ /www. stsci.edu/~postman/mdr_figures 



Fig. 6. — Image cutouts of 21 original (2 = 0.83) and redshifted (z = 1.24) MS1054-0321 galaxies and their corresponding morphological 
classifications. The original image is on the left. The galaxies shown here are spectroscopically confirmed cluster members. The redshifted 
images are constructed to match the exposure level used for our zs50 mosaic of RDCS1252-2927 (see TableQJ. 



ACS exposure of MS1358+6245 (z = 0.33) that were 
in common with the extensive study of this system per- 
formed by Fabricant et al. (2000). Agreement between 
the MP classifications and those from Fabricant et al. was 
achieved ~ 80% of the time with no systematic bias seen 
in the discrepant classifications (see Figure [7|). We thus 
conclude that our E, SO, and Sp+Irr classification scheme 
is robust and produces Hubble types that are compara- 
ble in accuracy to visual morphological data used in other 
studies. 

3.2. Field Morphological Population Fractions 

Figure El shows the E, SO, and Sp-flrr fractions as a 
function of zgso for galaxies either in the field. The data 
points are the fractions derived from our ACS data and 
our visual classifications. The results in this figure are 
based exclusively on galaxies with spectroscopic or pho- 
tometric redshifts that are incompatible with their be- 



ing cluster members. The grey bands show the typical 
range in low-density population fractions derived from 
local galaxy redshift surveys {e.g., PG84, D97, Goto et 
al. 2003a). The local {z < 0.1) and distant (0.5 < z < 1) 
field (i.e., low-density) galaxy populations appear to have 
similar fractions of E, SO, and Sp-|-Irr systems. 

4. LOCAL PROJECTED GALAXY DENSITY ESTIMATION 

We compute the local projected galaxy density in two 
different ways to ensure robustness - the nearest N neigh- 
bors approach used by D80 and D97 and a friends-of- 
friends algorithm. Both methods yield consistent results 
and we therefore present our results in terms of the near- 
est N neighbors based density unless otherwise specified. 
Appendix fXlprovides a demonstration of the consistency 
of these two density estimation techniques. In the near- 
est neighbor method, one computes the area of the region 
containing the N nearest neighbors and then derives the 
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Fig. 7. — Left Panels: The population fractions of E and SO galaxies as a function of 1775 magnitude for the cluster MS1358+6245 
at 2 = 0.33. The population fractions shown are from Fabricant et al. (2003) (based on their WFPC2 mosaic), and from this paper 
(based on the 1775 ACS WFC image and a redshifted version of this image out to 2 = 0.83). There are no obvious systematic offsets 
between our fractions and those in Fabricant et al. nor any systematic changes as we redshift the data out to 2 = 0.83. The dashed lines 
show the population fractions obtained by averaging the results of all classifiers. Right Panels: The population fractions in the original 
MSf054-0321 ACS image and in a version redshifted to 2 = f .24. The population fractions at each magnitude are offset from one another 
by a small amount for clarity. 



corresponding projected density at each galaxy location 
from the expression: 

„ fcorr{McL,MRef) f / n_1n \ 



, fc=i 



where Si is the projected galaxy density about a given 
galaxy, fcorriMcL, Mfj^f) is a correction factor that en- 
sures the density is always measured with respect to 
a common fiducial luminosity that corresponds to that 
used in low-redshift studies of the MDR, N is the num- 
ber of nearest neighbors, w(rnk, c^) is the selection func- 
tion (that can depend on magnitude and color - see Ap- 
pendix ^ for details), N^kgd is a background contamina- 
tion correction (if needed), ri^r is the solid angle of the 



region containing the N nearest neighbors (a rectangular 
region in our implementation), and Da is proportional to 
the angular diameter distance to the cluster (essentially 
the conversion between arcsec and projected Mpc). The 
correction factor is 



fcorr{McL , Mjief ) 



(2) 



where M^ef is the absolute magnitude limit to which 
we measure all densities, Mcl is the available limit 
for the cluster being analyzed, and $(M) is the galaxy 
luminosity function (with My{z = 0) = —21.28 and 
a — —1.22). We choose Mnef to match that of the 
original D80 study at z = {My = -19.27 for our 
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Full Resolution Figure Available at 

http:/ /www. stsci.edu/~postman/mdr_figures 



Fig. 8. — Color composite images of the 16 galaxies in common between our morphological sample for MS1054-0321 and those classified 
as Merger/Peculiar by van Dokkum et al. (2000). Our morphological classification of the object centered in each cutout is indicated in the 
lower left corner. The pictures here are made from the Veoei *775i a^nd ^850 ACS/WFC images. 



adopted cosmological parameters) but as we are sam- 
pling lookback times over which significant evolution in 
the characteristic magnitude of a Schechter luminosity 
function is detected, we allow Mjief to vary with red- 
shift as Mnefiz) = Miief{z = 0) — 0.8z. For our data, 
fcorr lies in the range [1.2,3.0]. We use the density de- 
rived from the N = 7 nearest neighbors but the results 
are not sensitive to this choice in the range 5 < iV < 10. 

The background correction is applied only when we are 
using samples requiring a statistical background subtrac- 
tion. For samples based on spectroscopic or photomet- 
ric redshifts, the background subtraction (if needed at 
all) is computed using the prescriptions described in Ap- 
pendix ^ Our statistical background correction, N^hgd, 
is derived from a combination of ACS and ground-based 
data. We use ACS observations of the HDF and Tadpole 
Galaxy (Benitez et al. 2004) to generate the surface den- 
sity of field galaxies when 1775 > 23.0. For 1775 < 23.0, 



we transform the number counts from the large /—band 
survey of Postman et al. (1998b) to the required ACS 
bandpasses. 

Density estimation is most accurate when using galaxy 
samples that are based on spectroscopic or photomet- 
ric redshift information as fore/background objects are 
effectively excluded from the analyses. In cases where 
sufficient redshift information is not available, one can 
use statistically-subtracted background corrected den- 
sity estimates. Such estimates are only reliable in dense 
regions (> 80 galaxies Mpc~^) where the cluster popula- 
tion dominates the counts. For reference, the la error in 
the field galaxy surface density at zs5q—24: is 5.4 galaxies 
arcmin"^, which corresponds to a projected density of 
26 and 22 galaxies Mpc"^ at z = 0.83 and z = 1.27, re- 
spectively. The statistically-subtracted background cor- 
rected population fractions and densities can be biased if 
there happens to be a significant overdensity in the line of 
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sight to the cluster. Thus, the most rehable resuhs are 
those based on samples with complete or nearly com- 
plete spectroscopic or photometric redshift information. 
In Appendix we demonstrate that a reliable compos- 
ite MDR or MRR can be derived from a combination 
of spectroscopic samples, photo-z samples, and samples 
with statistically subtracted background corrections pro- 
viding each such sample is limited to its optimal density 
regime. 

5. THE MORPHOLOGY - DENSITY AND MORPHOLOGY - 
RADIUS RELATIONS AT Z ~ 1 

We present our MDR and MRR results in Figures ITUl 
through!!^ Figurelinishows the composite MDR derived 
from all clusters in the sample using the best available 
data. The composite MDR is derived from the spectro- 
scopic samples for MS1054-0321 and RXJ0152-1357 for 



densities below 1000 galaxies Mpc ^ and their photo-z 
selected samples for densities above 1000 galaxies Mpc~^, 
the photo-z sample for RDCS1252-2927 (for aU densi- 
ties), and the statistically subtracted background results 
for CL1604-f4304, CL1604-f 4321, RDCS0910-f 5422, and 
RXJ0849-I-4452. The statistically subtracted back- 
ground results are used only when S > 80 galaxies 
Mpc~^ (roughly 3 times the amplitude of the typical fluc- 
tuations in the surface density of fore/background galax- 
ies; see section 0)). As in Sm04, we find that the MDR 
exists at z ~ 1. The z I MDR results from Sm04 
are shown for comparison. Our /e+so vs. local den- 
sity results are consistent with those of Sm04 to within 
the la uncertainties: we also find a less rapid increase in 
/e+so with increasing density than is seen at low redshift. 
However, the elliptical fraction, /e, shows no significant 
departures from the low redshift /E-density relation — 
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although, as demonstrated in Figure 01 the classification 
uncertainties in our /e and /so values are about twice as 
high as the uncertainties in our /e+so and /sp+irr values. 
The most notable difference between our current results 
and those at z < 0.2 (e.g., see Figure 9 in Fasano et al. 
2000) is the significantly lower fraction of SO galaxies: 
averaged over all densities with E > 30 galaxies Mpc~^, 
our mean /so= 0.20 ±0.12 where the error includes both 
the uncertainties from counting statistics (±0.035) and 
classification errors (±0.11). The typical low-z /so value, 
when averaged over the same density range, is 0.46 ±0.06 
(D80, D97, PG84) - a factor of - 2 higher than what is 
found at z ~ 1. A significant decline in the SO pop- 
ulation, relative to that seen in the current epoch, has 
previously been reported at redshifts as low as z ~ 0.4 
(D97, Fasano et al. 2000). Furthermore, over the range of 
densities being probed in this study, /so exhibits only a 
weak dependence on the projected density, analogous to 
what is seen at similar projected densities at lower red- 
shifts {e.g., D97). The shallower growth of /e+so with 
increasing density seen at z ~ 1 by Sm04 and by us thus 
appears to be due to a significant deficit of SO galax- 
ies and an excess of spiral galaxies relative to similar 
environments in the current epoch. This provides fur- 
ther support for observations suggesting that it is the SO 
and spiral population fractions that are experiencing the 
most significant changes with time over the past 8 Gyr 
{e.g., Moss & Whittle 2000; Fasano et al. 2000; Kodama 
& Small 2001; Treu et al. 2003; Sm04). 

Figure ^2 shows the morphology - radius relation 
(MRR) for the same composite cluster sample used in 
deriving the MDR in Figure ^| For each cluster we 
have computed an estimate of r2oo, the radius contain- 
ing a mean overdensity of 200 relative to a critical uni- 
verse, based on equation 8 in Carlberg, Yee, & EUingson 
(1997). The derived r2oo values are listed in Table [3 
For RDCS0910+5422 and RXJ0848+4452, we assume 
a = 750 km s~^. We still need to define a cluster cen- 
ter for each system and for that we use the centroid of 
the X-ray surface brightness distribution. X-ray imag- 
ing in which the hot ICM is detected with adequate S /N 
levels is available for all clusters except CL1604+4321. 
For CL1604+4321, we use the centroid of the distribu- 
tion of spectroscopically confirmed members. In the case 
of RXJ0152-1357, the X-ray distribution shows two well 
separated peaks (Delia Ceca et al. 2000) and we thus sub- 
divide the data for this cluster into two separate samples 
- a NE and SW component - and measure the radial dis- 
tance of cluster galaxies in each subsample relative to the 
nearest X-ray peak. The low-z reference for the MRR 
is taken from Whitmore & Gilmore (1991), with their 
radii converted to r2oo units using our cosmology and as- 
suming a mean cluster redshift of z — 0.04 and a mean 
cluster velocity dispersion of 750 km s^^. To aid in the 
comparison between the MDR and the MRR, we provide 
the approximate projected density - radius relationship 
log]^Q(S/635) ~ — 1.63(r/r2oo), which is derived from our 
data. This approximate relation is not particularly accu- 
rate for radu less than 0.2r2oo and the cluster-to-cluster 
variation about this relation is substantial (factors of 2—4 
variation in the projected density at a given r2oo-scaled 
radius), which is not surprising given the asymmetric 
galaxy distributions in many of our z > 0.8 clusters {e.g., 
Figure^ . The key features of the MRR at z ~ 1 are that 



(1) the bulk of the transition from a /sp+iir consistent 
with that in the field environment to its minimum value 
occurs within 0.6 r2oo (0.6 x r2oo corresponds to physical 
radii of 550 kpc - 1.1 Mpc for these clusters), (2) the 
z I Je+so value, at a given radius, is systematically 
less than the low-z /e+so for (r/r2oo) ^ 1-0, and (3) the 
/so - radius relation shows the most significant difference 
from the current epoch relationship. All of these charac- 
teristics are consistent with those inferred from the z ~ 1 
MDR. Given the significant asymmetry of the galaxy dis- 
tributions in some of these clusters, however, it is likely 
that the MRR is being diluted and is, thus, not as clean 
an indicator of morphological population gradients as the 
MDR for this particular cluster sample. 

5.1. Evolution of the Morphology-Density Relation 

An additional way to assess the presence of SO galaxies 
in clusters is to characterize the E+SO ellipticity distri- 
bution, as was done by D97. As shown earlier in Fig- 
ure |S1 the ellipticity distributions of lenticular and ellip- 
tical galaxies differ substantially. For the full ACS galaxy 
sample shown in Figure[Sl the hypothesis that the E and 
SO ellipticity distributions are drawn from the same par- 
ent population is rejected at greater than the 99.999% 
confidence level. Thus, if the z '--^ 1 bulge-dominated 
cluster galaxy population were truly devoid of a signifi- 
cant population of SO galaxies, the ellipticity distribution 
of the E+SO galaxies would more closely resemble that of 
pure ellipticals and would not include a significant com- 
ponent of objects with ellipticities beyond 0.5. Figure [T^ 
shows the ellipticity distributions of the E+SO cluster 
galaxy populations for three of our z > 0.8 cluster galax- 
ies along with galaxies from five clusters with redshifts 
in the range 0.25 < z < 0.55. To make a fair compari- 
son of these two cluster samples, we must select objects 
that lie in similarly dense environments. Therefore, the 
galaxies used in this comparison are all selected from en- 
vironments where the local projected density is > 100 
galaxies Mpc~^. The z > 0.8 sample used here con- 
sists of the photo-z selected cluster members in RXJ0152- 
1357, MS1054-0321, and RDCS1252-2927. We limit the 
z > 0.8 cluster sample to the three clusters with good 
photo-z data to ensure we are selecting probable cluster 
members. The z < 0.6 cluster data are from GTO ACS 
observations of 5 strongly lensing clusters in the range 
0.25 < z < 0.55 (HST FID #9292). The 5 clusters are 
Zwl455+2232 (z = 0.258), MS1008-1224 (z = 0.301), 
MS1358+6245 (z = 0.328), CL0016+1654 (z = 0.54), 
and MS J0454-0300 (z ^ 0.55). Visual classification of 
all galaxies with 1775 < 22.5 in each z < 0.6 cluster ACS 
image was performed by MP. The rest-frame wavelengths 
sampled by the 1775 filter for 0.25 < z < 0.55 typically lie 
in the 1^— band, thus any "morphological" k-corrections 
between this sample and the z > 0.8 sample should be 
small. A total of 798 galaxies were classified in the five 
z < 0.6 clusters. As the low-z ACS images were all sin- 
gle pointings centered on the cluster core, the vast ma- 
jority of the E and SO galaxies identified are likely to 
be cluster members. For reference, the ACS/WFC sub- 
tends 800 kpc at z = 0.26 and 1.3 Mpc at z = 0.55. 
A Kolmogorov-Smirnov test finds that the distribution 
of ellipticities of E+SO galaxies in the z > 0.8 clusters 
is inconsistent with being drawn from a pure elliptical 
ellipticity distribution at the 97.0% confidence level. A 
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Fig. 10. — The MDR for all clusters based on the best available data. See text for details. Errorbars include the uncertainties in counting 
statistics and morphological classification. The previous results from low-2 surveys (D80, D97, PG84, SDSS) and from 2 ~ 1 (Sm04) are 
shown for reference. The triangular markers along the y-axis show the low-density population fractions from the SDSS (Goto et al. 2003a). 



Wilcoxon rank-sum test, which is better suited to mea- 
suring differences in the mean vaiues of two distributions 
than a KS test, finds the E-fSO and E gaiaxy distribu- 
tions for the z > 0.8 ciusters differ at the S.lcr ievei. We 
thus reject the hypothesis that there are no SO galaxies in 
dense environments at z ~ 1. The E-I-SO and E galaxy 
ellipticity distributions for the 0.25 < z < 0.55 cluster 
sample are inconsistent with each other at the 99.998% 
confidence level. 

We can also apply a robust test {i.e., one that primar- 
ily relies on our ability to distinguish only between E-I-SO 
and Sp-|-Irr) in an attempt to constrain the evolution of 
the SO population by comparing the z > 0.8 E-I-SO el- 
lipticity distribution with that for the 0.25 < z < 0.55 
cluster galaxy sample. The cumulative distribution func- 
tions for the E-I-SO ellipticities in the two different cluster 
samples are shown in the top panel in Figure The 
ellipticity distributions for the 0.25 < z < 0.55 and the 



z > 0.8 cluster samples are, however, completely con- 
sistent with one another (KS test rejects inconsistency 
only at the ~ 10% confidence level), suggesting that, at 
least in these particular samples, any evolution in the 
relative E and SO population fractions does not manifest 
itself as a significant difference in the E-I-SO ellipticity 
distribution. It is possible that with larger and more 
homogeneously selected cluster samples over a range of 
redshifts this test would provide a quite objective way to 
measure the evolution of the morphological population 
fraction in dense environments. 

The Je+so as a function of lookback time and local 
density is shown in Figure [T51 This figure is modeled af- 
ter Figure 3 in Sm04. The Sm04 results are reproduced 
for reference in the lower panel of Figure 1131 The upper 
panel of Figure El shows the analogous results for our 
ACS z > 0.8 and 0.25 < z < 0.55 cluster samples. The 
results in the upper panel are based on the mean popu- 
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Fig. 11. — The Morphology — Radius relation for all clusters based on the best available data. See text for details. Errorbars include the 
uncertainties in counting statistics and morphological classification. Cluster centers are determined from the centroid of the X-ray surface 
brightness distribution except for CL1604+4321, where the centroid of the distribution of confirmed spectroscopic members is used instead. 
The solid line is the low-2 morphology — radius relation from Whitmore & Gilmore (1991), converted to r200 units assuming a mean cluster 
redshift of z = 0.04 and a mean cluster velocity dispersion of 750 km s~^. The triangular markers along the y-axis show the low-density 
population fractions from the SDSS (Goto et al. 2003a). 



lation fractions within logarithmic density bins 0.4 units 
wide centered at 10, 100, and 1000 h^^ galaxies Mpc^'^. 
Our /e+so values are in good agreement with the Sm04 
values - any differences are comparable with the uncer- 
tainties. We corroborate the key observational result of 
the Sm04 study that the most significant differences be- 
tween the MDR at low redshift and high redshift are 
confined to regions where the projected galaxy density is 
larger than ~ 40 galaxies Mpc^^. While the agreement 
is re-assuring given the significant overlap of the clus- 
ters used in the two programs, and bolsters the concept 
of using visual morphological classification in compara- 
tive studies at high-z, there is also significant room for 
reducing the existing uncertainties that will be achieved 
only when a much (~ 10 x) larger sample of z ~ 1 cluster 



galaxy morphologies is available. Our derived population 
fractions as a functions of projected density and radius 
are also provided in Tables O and 0] for easy reference. As 
noted above, our results suggest that the observed MDR 
evolution is primarily driven by evolution in the fractions 
of SO and Sp-|-Irr galaxies. 

5.2. A Correlation between /e+so Oind ? 

Our sample exhibits a potentially interesting trend be- 
tween the early type population fraction and the clus- 
ter bolometric X-ray luminosity. The bulge-dominated 
galaxy population fraction in clusters with high X-ray 
luminosity (Lx,Boi> 1-5 x lO''^^ h^Q erg s~^) is higher 
than that in low X-ray luminosity (ix.Boi^ 3 x 10^'^ 
erg s~^) clusters, even at the highest projected densities. 
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Fig. 12. — The distribution of ellipticities for the bulge-dominated galaxies that are likely cluster members in 2 different redshift ranges. 
The lower redshift cluster data are also based on ACS imaging and come from observations of 5 strongly lensing clusters in the range 
0.25 < z < 0.55. See section for details. The z > 0.8 galaxies are from our photometric redshift sample for the clusters RXJ0152- 
1357, MS1054-0321, and RDCS1252-2927. The E+SO ellipticity histograms are drawn with thick hues. The elliptical galaxy ellipticity 
distributions are represented by the light grey shaded histograms. The uppermost panel shows the cumulative distribution functions for 
the ellipticities of the E+SO galaxies in each cluster sample. 



This is demonstrated in Figure [HI Here we show /e+so, 
/e, and /so as a function of the bolometric X-ray lumi- 
nosity for the galaxies within r2oo for each cluster. D80 
found a modest increase in /so, a corresponding decrease 
in /sp+irr, and no change in /e in his sub-sample of 8 
high X-ray emitting clusters relative to the full sample 
of 55 clusters. Least square, error- weighted, fits to the 
trends in Figure IT4l give 



/ 



E+SOCX-^x,Bol 



0.33±0.09 



/eocLx,Bo1 



0.15±0.09 



0.18±0.09 



/soOf-^^x,Bol 

The linear correlation coefficients for the fits to the 
log(Lx,Boi)- /e+so, log(Lx,Boi)- /e, and log(Lx3oi 



)- /so data are 0.82 (97.6% CL), 0.84 (98.2% CL), and 
0.75 (94.8% CL), respectively. The correlations are sig- 
nificant at the 2 — 3a level. 

A correlation between /e+so and could, in prin- 
ciple, be produced either as a consequence of environ- 
mental interactions subsequent to the formation of the 
cluster galaxies and/or as a result of the initial condi- 
tions present at the time of their formation. In a sim- 
plified example of the former scenario, the ram-pressure. 
Pram, acting on a galaxy moving through the intracluster 
medium (ICM) is proportional to pv^ . If the ICM is in 
hydrostatic equilibrium, then the bolometric X-ray lumi- 
nosity, Lx, is proportional to T^^'^ B? {e.g., Ettori et 
al. 2004). The total mass of the system, yitot, is propor- 
tional to Ty^^^ and is also proportional to v^. The same 



The Morphology-Density Relation in 1 Clusters 



17 



0.8 - 



" 0.6 



0.4 



0.2 
1 



0.8 



03 

J> 0.6 



0.4 - 



0.2 



This Work 



»S<10 Mpc"' 
>S=100 Mpc"^ 
I E=1000 Mpc"' 



Smith et al. 2004 



I — 



_ - 



1 

1 



^ 



-\ ^K^- ^ 



_ _ J _ - 



Lookback 



6 

(Gyr) 



10 



Fig. 13. — Bottom: The evolution of the MDR reported by Sm04. Loolcbaclc times liavc been transfornied to compensate for the small 
difference between the cosmology in Sm04 and tiiat adopted in this paper - specifically, our clioice of a slightly larger Hubble constant (70 
vs. 65). The dashed lines representing the change in /e+so with time at the 3 different density regimes are reproduced in all three panels 
for reference. Top: Our derived E+SO fractions at projected densities corresponding to < 10, 100, and lOOOhgg galaxies Mpc~^ (to match 
Sm04) for our combined ~ 1 cluster sample and our 0.25 < z < 0.55 cluster sample. The latter sample is not suitable for measuring 
low-density population fractions. 



family of scaling laws gives Ly^(xT^ ^ although observa- 
tionally a steeper relation is found where L^(xTJ^'^ {e.g., 
Ponman et al. 1996; Mulchaey & Zabludoff (1998); Xue 
and Wu (2000)). These relations can be manipulated 
to yield Pram C)cLx^ °^°'^! where the exponent value de- 
pends on the exponent in the Lx — relation. The ram- 
pressure is stronger within more luminous X-ray clusters. 
Therefore, if ram-pressure stripping were the dominant 
mechanism responsible for the origin of SO galaxies in 
clusters, /e+so should exhibit a positive correlation with 
Lx- However, many previous observations, including the 
relatively weak dependence of /so on projected density, 
strongly suggest that ram pressure stripping alone cannot 
explain the morphological mix in clusters {e.g., D80; D97; 
Kodama & Small 2001; Okamoto & Nagashima 2003). 



Alternatively, a positive correlation between /e+so and 
Lx could be produced as a consequence of simple dynam- 
ics - the most massive (and, hence, most X-ray luminous) 
clusters will collapse earlier and any environmentally- 
driven processes that produce bulge-dominated galaxies 
will therefore have been active for a longer period of time 
at any subsequent redshift. This simple timescale argu- 
ment could result in a positive correlation between the 
fraction of early-type galaxies and observational proxies 
for cluster mass, regardless of whether or not the effi- 
ciencies of the transformation processes are correlated 
with the properties of the ICM. However, this would only 
be true if the timescales required to establish a signifi- 
cant population of early-type galaxies were long com- 
pared to the collapse time. If the timescale for estab- 
lishing the early-type population were comparable to the 
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collapse time (i.e., if deeper potential wells are "born" 
with a higher early-type population fraction), then the 
/e+so — relation may be telling us more about cluster 
and galaxy formation processes than about the cluster 
evolution process. 

As the results in Figure El are based on only 7 
clusters and as the significance of the correlations are 
only significant at the < 3(t level, we refrain from over- 
interpretation. Indeed, the correlations between /e+so 
and Tx and /e+so and cluster velocity dispersion, cr, are 
not significant (see Figure I15|l . In the literature, the 
mass dependence of cluster galaxy evolution has been 
a subject of debate. For example, Fairley et al. (2002) 
studied eight X-ray selected clusters and did not find any 
dependence of the blue galaxy fraction on Lx - a trend 
that might be expected if /e+so correlated with Lx- De 
Propris et al. (2004), using a larger sample of 60 clus- 
ters from the 2dF Galaxy Redshift Survey, finds that the 
blue galaxy fraction does not depend on the velocity dis- 
persion of the cluster galaxies. However, there is tenta- 
tive evidence that cluster integrated star formation rates 
{e.g., Finn et al. 2004) correlate with Tx and (Homeier 
et al. 2005), which is presumably a dependence on cluster 
mass. Furthermore, Zabludoff & Mulchaey (1998) found 
that groups of galaxies exhibit a strong correlation be- 
tween the bulge-dominated galaxy fraction, /e+soj and 
the group velocity dispersion. And while the specific re- 
lation they found cannot be extended to very massive 
clusters (their predicted /e+so reaches a value of unity 
at a velocity dispersion of ~ 700 km s~^), it does sug- 
gest that groups of galaxies have a positive correlation 
between /e+so and Lx through the observational relation 
Lxoc (Mulchaey & Zabludoff 1998). Margoniner et 
al. (2001) and Goto et al. (2003b) find that blue fractions 
of cluster galaxies are lower for richer clusters, a result 
that is consistent with a positive /e+so ~ relation- 
ship. 

The significance of the /e+so ~ relation clearly 
needs to be studied using a larger sample so that sub- 
division by mass (or a suitable observational proxy for 
mass) can be conducted for a far greater number of clus- 
ters. The HST/ACS snapshot program of 73 homoge- 
neously selected X-ray clusters (ID#10152, M. Donahue, 
P.I.) should provide the sample needed to assess this 
relation in the range 0.3 < z < 0.7. If a significant 
/e+so ~ relation is ultimately found, one implication 
is that the MDR and its evolution may be dependent on 
the total cluster mass. Certainly on galaxy-mass scales, 
there is evidence that the star formation process occurs 
more rapidly in systems with higher stellar mass {e.g., 
Heavens et al. 2004). For the present sample, correla- 
tions between the fraction of early type galaxies and the 
X-ray properties of the clusters provide, at best, a strong 
inspiration to study the trends with samples explicitly 
geared to investigate this relationship. 

6. DISCUSSION 

As originally discussed by TuUy and Shaya (1984), the 
measurement of the evolution in the MDR is a significant 
step in understanding the relative roles of environment 
and initial conditions in establishing morphological pop- 
ulation gradients. The less rapid growth in the bulge- 
dominated galaxy population with increasing density (or 
decreasing radius) that is seen at 2: ~ 1, implies that en- 



vironment must play an important role in the establish- 
ment of the current epoch MDR. Sm04 explore a range 
of morphological transformation scenarios based on the 
expression: 

/S0,z=1-/E+S0,z=i - /E,2=o.5 -ITf ^ — [6) 

where /so z=i the SO population fraction in clusters 
at z = 1, /e+so z=i E-I-SO fraction in clusters at 

z = 1, /e z=o 5 the elliptical fraction in clusters at 
z — 0.5, Nz=Q.5/Nz=i is the ratio of the number of galax- 
ies in clusters at z = 0.5 to that at z = 1, and AA'^ is 
the change in the number of cluster ellipticals between 
z — I and z = 0.5. By combining varying amounts of in- 
fall (including infall of early-type galaxies), mergers, and 
cannibalism, Sm04 find that for reasonable levels of each 
process the /so at z = 1 is typically less than 0.1. The 
exception is for the case where there is no infall and 10% 
of the spiral galaxies merge in pair-wise manner to form 
ellipticals. In this case, the model predicts /so= 0.18 at 
z = 1. The degree to which our observations conflict 
with the predictions for very low SO population fractions 
in z = 1 clusters is, of course, directly tied to the uncer- 
tainty in our /go measurements. Our mean SO population 
fraction is inconsistent with being zero at about the 90% 
confidence level, which is not strongly in conflict with the 
Sm04 predictions for /so z=i- However, given that the el- 
lipticity distribution of the E-I-SO galaxies in the z ~ 1 
clusters is inconsistent with being drawn from a popu- 
lation dominated largely by pure elliptical morphologies 
and given that we do get a mean /so of 0.20 ± 0.12 for 
S > 30 galaxies Mpc~^, we investigate minor variations 
to the Sm04 scenario that can increase in the predicted 
SO population fraction at z ~ 1. An obvious choice is to 
explore reasonable modifications of the value of /e z=o 5- 
Sm04 adopted /e,z=o.5= 0.6. However, we find that for 
S > 100 galaxies Mpc~^ a value in the range 0.4 — 0.5 ap- 
pears to be closer to what the data suggest. Decreasing 
/e,z=o.5 to 0.5 increases the range of the Sm04 predic- 
tions to 0.09 </so z=i— 0-28 for the models considered. 
This would make their predictions and our mean /so val- 
ues consistent within la. However, doing so does change 
the conclusions of Sm04 a bit - namely that ~50% of the 
lenticulars in clusters could already be in place at z = 1 
and that the remaining half form largely between z ^ 0.5 
and the present epoch. 

Different morphological transformation mechanisms 
operate in different environments and, thus, by identify- 
ing the radius or density where the morphology of cluster 
galaxies start to change, we can hope to identify the un- 
derlying physical mechanisms. The breaks in the MDR, 
first characterized by PG84 and more recently by Goto 
et al. (2003), suggest different processes are probably re- 
sponsible for the origin of elliptical and SO galaxies. For 
example, ram-pressure stripping (Gunn & Gott 1972) is 
only efficient within the cluster core (typically < 250 
kpc) where the ICM density is high enough to allow the 
dynamic pressure to overpower the gravitational restor- 
ing force of interstellar gas in a galaxy's disk. Figure ITT] 
indicates that the fraction of elliptical galaxies starts to 
increase inwards of 0.6 r2oo, or in terms of local galaxy 
density, at around 70 galaxies Mpc~^ (Figure EU- This 
environment matches that where the population frac- 
tion of Sp-|-Irr galaxies starts to decline. The scale of 
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Fig. 14. — The dependence of the bulge-dominated galaxy population fractions, /e+so, /ej and /so, on the bolometric cluster X-ray 
luminosity within a radius corresponding to r2oo- Names of the individual clusters are shown in the upper panel. The solid line is the best 
fit relation when the data are weighted by the inverse square of their uncertainties. The dashed line is the best fit when each data point is 
given equal weight. The dotted line is the best error-weighted fit with CL1604+4321 excluded. The data for RXJ0152-1357 in this figure 
includes both the NE and SW components of the cluster. Errorbars include the uncertainties in counting statistics and morphological 
classification. 



0.6 r2oo for the onset of the transition between the low- 
density field population and that in high density regions 
is in agreement with previous low redshift observations 
{e.g., Kodama et al. 2001; Treu et al. 2003; Goto et al. 
2003b). At distances beyond 250 kpc or so, ram pressure 
becomes much less effective as a morphological transfor- 
mation process and thus other environmental processes 
must play a significant role in the establishment of the 
MDR and MRR {e.g., Treu et al. 2003). Galaxy-galaxy 
merging, for example, is most effective when relative ve- 
locity of galaxies are comparable to the rotation of galaxy 
itself (~ 200 km s^^) and is thus likely to be most ef- 
fective in regions far from the cluster core. Galaxies be- 
yond 5 Mpc from the cluster center will, at ^; --^ 1, not 
yet have had sufficient time to traverse the cluster core 
and should, therefore, not experience tidal stripping or 



tidal triggering of the star formation due to the cluster 
gravitational potential. 

The lack of any significant redshift dependence in the 
/e - density relation coupled with the significant changes 
with redshift seen in the /so - density relation suggest 
that the origin of the MDR and MRR in high density 
regions is most likely determined by two distinct pro- 
cesses: (i) formation of elliptical galaxies at the clus- 
ter core, which occurred at redshifts significantly greater 
than unity and (ii) the formation of SO galaxies, which 
appears, for at least half of the lenticulars, to be a process 
that is still underway even at redshifts as low as z ~ 0.5. 
Evidence for different formation timescales for elliptical 
and SO galaxies comes from other types of studies as 
well. For example, the analysis of the color-magnitude 
relation (OMR) in RDCS1252-2927 by Blakeslee et al. 
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Fig. 15. — The bulge-dominated galaxy population fraction, /b+sOi within a radius corresponding to r2oo as a function of the X-ray gas 
temperature and cluster velocity dispersion. Names of the individual clusters are shown. Errorbars include the uncertainties in counting 
statistics and morphological classification. 



(2003) finds that the observed scatter in the CMR for 
SO galaxies in that system is consistent with relatively 
recent or even on-going star formation activity whereas 
the CMR for elliptical galaxies in this z = 1.24 cluster is 
compatible with no significant star formation activity for 
more than a Gyr. In addition, Poggianti et al. (1999) and 
Goto, Yagi, Tanaka, & Okamura (2004) derive shorter 
evolutionary timescales for blue to red spiral transfor- 
mation than for red spiral to red elliptical transforma- 
tion, consistent with the hypothesis that transformation 
of a galaxy's spectral energy distribution typically occurs 
more rapidly than morphological transformation. 

A fundamental question that remains is whether or 
not the MDR and MRR depend on luminosity. It is well 
established that intrinsically luminous, red galaxies are 
almost always elliptical and, hence, at some level limiting 
the analyses of the MDR to intrinsically luminous galax- 
ies may produce a different relationship between mor- 



phology and projected density than found in a sample 
containing a broader range of luminosities. Tanaka et al. 
(2005) demonstrate that there is indeed a change in the 
MDR at 2; < 0.07 as a function of luminosity based on 
an analysis of ^ 20, 000 SDSS galaxies. They find that 
galaxies fainter than M* -|- 1 do appear to have a dif- 
ferent MDR at low densities than brighter galaxies. We 
have divided our spectroscopic samples for MS1054-0321 
and RXJ0152-1357 into four apparent magnitude ranges 
^775 < 22, ^775 < 22.5, ^775 < 23, and 1775 < 23.5 corre- 
sponding to M* < M7*75 - 0.3, M* < M7*75 + 0.2, M* < 
M7*75 + 0.7, and M* < M^^^ + 1.2. We find no sig- 
nificant difference between the MDR derived for these 
4 samples. If the luminosity dependence of the MDR 
does not exhibit strong evolution, then our current sam- 
ple does not go deep enough at z = 0.83 to accurately 
measure the dependence. Our failure to detect any lu- 
minosity dependence in the MDR in our z = 0.83 spec- 
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troscopic sample is thus not inconsistent with the work 
done at low rcdshifts. Furthermore, our survey is not well 
suited to studying densities below 30 galaxies Mpc~^, 
well above the density where the Tanaka et al. (2005) 
find the strongest luminosity-dependent efii'ects. All we 
can conclude for now is that for S > 30 galaxies Mpc~^ 
and for M* < Aff^g + 1.2 the z = 0.83 MDR is not 
sensitive to luminosity. The caveat to this conclusion is 
that our ^775 data at z = 0.83 samples the rest-frame 
B— band. It may be that galaxies selected according to a 
redder rest-frame luminosity will show a diff'erent trend. 
The underlying goal of studying the MDR as a func- 
tion of luminosity is, of course, to assess the dependence 
of the MDR on galaxy mass. We are assembling NIR 
photometry for our cluster galaxies to provide more ac- 
curate stellar mass estimates and we will explore trends 
between morphology, density, r2oo, and galaxy mass in a 
future paper. 

7. CONCLUSIONS 

We have performed deep, multiband observations with 
the ACS/WFC of 7 clusters with 0.83 < z < 1.27 to 
study the morphological composition of the galaxy pop- 
ulation over 3 decades of local density and out to radii 
of up to 2 Mpc from the peak of the X-ray emission from 
the ICM. The key results are: 

1. The high sensitivity and angular sampling of the 
ACS/WFC enable reliable visual distinctions to 
be made between the major morphological galaxy 
classes (E, SO, Sp-hirr) down to ~ 0.25L/L* in the 
range 0.8 < z < 1.3. 

2. We confirm the results of Smith et al. (2004) that 
a morphology - density relation exists at 2; ~ 1. 
We also explicitly confirm the less rapid growth in 
/e+so with increasing density and that the change 
in the slope of the MDR is due primarily to changes 
in the high-density population fraction. A flatten- 
ing of the Je+so ~ density relation with increasing 
redshift can be a consequence of environmentally- 
driven transformation of galaxies from late to early- 
types. Our new results provide direct measure- 
ments of the lenticular population fraction at z ~ 1 
and we conclude that the observed difii'erences in 
the morphology - density relation at z ^ 1 are due 
primarily to a deficit of SO galaxies and an excess 
of Sp-|-Irr galaxies relative to the local galaxy pop- 
ulation. The /e - density relation does not appear 
to evolve over the range < z < 1.3. A deeper 
understanding of the implications of these results 
for models of galaxy and cluster formation will re- 
quire further exploration of the dependence of mor- 
phological population fractions on galaxy mass and 
an assessment of the frequency of cluster galaxy 
merger activity. 

3. The MDR at z = 0.83 is not sensitive to the rest- 
frame i?— band luminosity for galaxies with lumi- 
nosities brighter than M* + 1.2 and in regions with 
S > 30 galaxies Mpc~^. Work done at low redshift 
suggests, however, that luminosity effects may be 
more pronounced at fainter luminosities and in re- 
gions of lower density. Our present survey is not 
well-suited to studying these regimes. 



4. We directly measure the lenticular population frac- 
tion and find /so= 0.20 ± 0.12 when we average 
over densities with E > 30 galaxies Mpc~^. The 
error in /so includes contributions from both count- 
ing statistics (±0.035) and errors in our ability to 
visually classify lenticular galaxies (±0.11). Our 
z ^ 1 /so value is about a factor of 2 less than 
the /so seen in similarly dense environments in the 
local universe but is comparable to what is seen 
at 0.4 < z < 0.5. Our 20% population fraction of 
SO galaxies is higher than almost all of the scenar- 
ios proposed by Smith et al. (2004) to explain the 
shallower /e+so ~ density relation - they predict 
/so< 0.1. However, a small reduction in the ellipti- 
cal population fraction at z = 0.5 adopted by Smith 
et al. (2004) from /e= 0.6 to /e= 0.5, a value that 
appears to be in better agreement with the observa- 
tions, is sufficient to increase their predicted z = 1 
SO population fractions to 0.2 ± 0.1, overlapping 
our measurement. The distribution of ellipticities 
in z > 0.8 bulge-dominated (E+SO) cluster galaxies 
is inconsistent with an ellipticity distribution that 
would arise from a sample consisting solely of ellip- 
tical galaxies. In other words, our results suggest 
that rich clusters at z ~ 1 most likely have a signif- 
icant population of lenticular galaxies and, there- 
fore, a significant percentage of lenticular galaxies 
could have formed at redshifts z > 1.3. 

5. We measure the morphology - radius relation and 
find its evolution is consistent with that seen in the 
morphology - density relation: 

(a) the bulk of the transition from a /sp-i-irr con- 
sistent with that in the field environment to 

its minimum value occurs within a radius of 
0.6r2oo (which corresponds, on average, to 
densities > 70 galaxies Mpc~^ and to phys- 
ical scales less than ~ 750 kpc), 

(b) the z ~ 1 /e+so value, at a given radius, is 
systematically less than the low-z /e+so for 
radii less than ~ r2oo, and 

(c) the /so - radius relation shows the most sig- 
nificant difference from the current epoch re- 
lationship. 

However, elongation of and dumpiness in the 
galaxy distributions for many of our z > 0.8 
clusters makes interpretation of the azimuthally- 
averaged morphology - radius relation more dif- 
ficult (and perhaps less meaningful) than trends 
between morphology and local density. 

6. We find that the bulge-dominated galaxy popula- 
tion fraction, /e+sO) is mildly correlated with the 
bolometric X-ray luminosity of the cluster. Clus- 
ters with high X-ray luminosities have higher /e+so 
values within r2oo than clusters with lower X-ray 
Imninosities. In the present sample, the trend is 
significant at the < 3(7 level . A correlation be- 
tween /e+so and bolometric X-ray luminosity can 
arise as a consequence either of environmentally- 
driven transformation processes or initial condi- 
tions. However, we do not find significant corre- 
lations between /e+so and the X-ray temperature 
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or cluster velocity dispersion. A definitive study 
of the relation between galaxy population fraction 
and cluster X-ray properties will require the analy- 
sis of large homogeneously selected cluster samples. 
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TABLE 1 

Summary of Spectroscopic and HST Observations 



Cluster 


Rcdshift 




Nz.Ci, ACS 


Mosaic Area 
(sq.arcmin) 


Filter (Exp. Time in ksee) 


HST ID 


MS 1054-0321 


0.831 


327 


143 


35.5 


V6O6 (2.0), 4775 (4.0), Z850 


(4.0) 


9290, 9919 


RXJ 0152-1357 


0.837 


123 


93 


36.5 


'•625 (4.8), 4775 (4.8), ZgsO 


(4.8) 


9290 


CL 1604+4304 


0.900 


107 


20 


12.0 


Veoe (4.8), Isi4 (4.8) 




9919 


CL 1604+4321 


0.921 


130 


31 


12.0 


Veoe (4.8), hu (4.8) 




9919 


RDCS 0910+5422 


1.101 


~ 10 


~ 10 


12.2 


J775 (6.8), Z850 (11.4) 




9919 


RDCS 1252-2927 


1.235 


180 


31 


32.7 


J775 ( 7.2), 0860 (12.0) 




9290 


RXJ 0849+4452 


1.266 


90 


16 


33.7 


i775 (7.3), 0850 (12.2) 




9919 



TABLE 2 

Summary of Cluster X-ray and Kinematic Data 









Tx 


cr 


1'200 


X-ray Data 


Vel. Disp. 


Cluster 


(10" h-^ erg s- 


-^) (10" h^f erg s- 


1) (keV) 


(km g-^) 


(h-^ Mpc) 


Refcrcnee 


Reference 


MS 1054-0321 


7.78 ± 0.4 


16.43 ± 0.8 




1153+80 


1.79 


1, 2 


10 


RXJ 0152-1357 


5.74 ± 0.6 


18.57 ± 1.9 




919+168 


1.42 


2, 3 


11 


CL 1604+4304 


0.86 ± 0.13 


2.01 ± 0.3 


9 ci+l.bs 
■'•°^-0.69 


989+?^ 


1.48 


4 


12 


CL 1604+4321 


< 0.7 






649lfe 


0.95 




12 


RDCS 0910+5422 


0.78 ± 0.09 


2.14 ± 0.3 


7.20+11 




0.99° 


5 




RDCS 1252-2927 


1.90 ± 0.3 


6.60 ± 1.1 


6.50 ± 0.5 




0.94 


6, 7 


13 


RXJ 0849+4452 


1.41 ± 0.3 


2.85 ± 0.6 


5.80t5;^ 




0.91" 


8, 9 





^i'200 based on assumed cr = 750 km s 

^Gioia et al. 2004; Donahue 2004, priv. comm. 

2Romcr et al. 2000 

^Della Ceca et al. 2000 

"'Lubin et al. 2004 

^Stanford et al. 2002 

'^Rosati et al. 2004 

■^Lombardi et al. 2005 

SRosati et al. 1999 

^Stanford et al. 2001 
lOGioia et al. 2004 
llDcmaico et al. 2004 
l^Gal & Lubin 2004 
i^Demarco et al. 2004a 
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TABLE 3 

Population Fractions as a Function of Projected Density 



Cluster or 




logio(S) 










Sample 


N(S) 


(Gals Mpc"^) 


/e+so 


/e 


/so 


/Sp+Irr 


z ~ 1 composite 


136 


1.30 


0.45 ± 0.11 


0.32 ± 0.14 


0.13 ± 0.14 


0.56 ± 0.11 




220 


1.60 


0.35 ± 0.09 


0.22 ± 0.13 


0.13 ± 0.13 


0.65 ± 0.09 




237 


1.90 


0.53 ± 0.09 


0.31 ± 0.13 


0.22 ± 0.13 


0.47 ± 0.09 




138 


2.20 


0.62 ± 0.10 


0.39 ± 0.14 


0.23 ± 0.14 


0.38 ± 0.10 




87 


2.50 


0.71 ± 0.12 


0.47 ± 0.15 


0.24 ± 0.15 


0.29 ± 0.12 




24 


2.80 


0.65 ± 0.21 


0.53 ± 0.23 


0.12 ± 0.23 


0.35 ± 0.21 




14 


3.10 


0.87 ± 0.27 


0.71 ± 0.29 


0.16 ± 0.29 


0.12 ± 0.27 


z ^ 1 composite 


185 


1.06 ± 0.2 


0.30 ± 0.10 


0.24 ± 0.13 


0.06 ± 0.13 


0.70 ± 0.10 




595 


2.06 ± 0.2 


0.51 ± 0.07 


0.31 ± 0.12 


0.20 ± 0.12 


0.49 ± 0.07 




38 


3.06 ± 0.2 


0.73 ± 0.17 


0.59 ± 0.20 


0.14 ± 0.20 


0.27 ± 0.17 


z < 0.6 composite 


435 


2.06 ± 0.2 


0.63 ± 0.08 


0.34 ± 0.12 


0.29 ± 0.12 


0.37 ± 0.08 




66 


3.06 ± 0.2 


0.80 ± 0.14 


0.50 ± 0.17 


0.30 ± 0.17 


0.20 ± 0.14 


MS 1054-0321 


79 


> 2.00 


0.80 ± 0.13 


0.49 ± 0.16 


0.31 ± 0.16 


0.20 ± 0.13 


RXJ 0152-1357 


63 


> 2.00 


0.81 ± 0.14 


0.56 ± 0.17 


0.25 ± 0.17 


0.20 ± 0.14 


CL 1604+4304 


62 


> 2.00 


0.53 ± 0.14 


0.33 ± 0.17 


0.20 ± 0.17 


0.47 ± 0.14 


CL 1604+4321 


86 


> 2.00 


0.39 ± 0.12 


0.25 ± 0.15 


0.14 ± 0.15 


0.61 ± 0.12 


RDCS 0910+5422 


43 


> 2.00 


0.50 ± 0.16 


0.37 ± 0.19 


0.13 ± 0.19 


0.50 ± 0.16 


RDCS 1252-2027 


k; 


> 2.(10 


(1.8.5 1 O.Ki 


0.62 + n.l8 


0.23 + 0.18 


0.15 + 0.16 


RXJ U>M!)+ll."i2 


lUl 




().:!.-) 1 U.ll 


0,:',:! 1 0.11 


0.02 ± O.il 


0.05 ± 0.11 



Errors are the quadrature sum of counting statistics and classi- 
fication uncertainty. See text for details. 
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TABLE 4 

Population Fractions as a Function of R200 Radius 



Cluster or 
Sample 


N(r/r2oo) 


r/r200 


/e 


+ S0 


/e 


/so 


/Sp + I 




2 ~ 1 composite 


63 


0.05 


0.75 


± 


0.14 


0.52 


± 


0.17 


0.23 


± 


0.17 


0.26 


± 


0.14 




80 


0.15 


0.52 


± 


0.13 


0.37 


± 


0.16 


0.15 


± 


0.16 


0.48 


± 


0.13 




90 


0.25 


0.56 


± 


0.12 


0.35 


± 


0.15 


0.21 


± 


0.15 


0.44 


± 


0.12 




83 


0.35 


0.45 


± 


0.13 


0.24 


± 


0.16 


0.21 


± 


0.16 


0.55 


± 


0.13 




94 


0.45 


0.42 


± 


0.12 


0.25 


± 


0.15 


0.17 


± 


0.15 


0.58 


± 


0.12 




99 


0.55 


0.27 


± 


0.12 


0.21 


± 


0.15 


0.06 


± 


0.15 


0.73 


± 


0.12 




114 


0.65 


0.32 


± 


0.11 


0.20 


± 


0.15 


0.12 


± 


0.15 


0.69 


± 


0.11 




129 


0.75 


0.30 


± 


0.11 


0.15 


± 


0.14 


0.15 


± 


0.14 


0.70 


± 


0.11 




96 


0.85 


0.35 


± 


0.12 


0.11 


± 


0.15 


0.24 


± 


0.15 


0.64 


± 


0.12 




74 


0.95 


0.25 


± 


0.13 


0.15 


± 


0.16 


0.10 


± 


0.16 


0.75 


± 


0.13 




23 


1.05 


0.32 


± 


0.22 


0.32 


± 


0.24 


0.00 


± 


0.24 


0.68 


± 


0.22 




11 


1.15 


0.19 


± 


0.31 


0.00 


± 


0.32 


0.19 


± 


0.32 


0.81 


± 


0.31 


MS 1054-0321 


130 


< 1.0 


0.78 


± 


0.11 


0.46 


± 


0.14 


0.32 


± 


0.14 


0.22 


± 


0.11 


RXJ 0152-1357 


125 


< 1.0 


0.72 


± 


0.11 


0.46 


± 


0.14 


0.26 


± 


0.14 


0.27 


± 


0.11 


CL 1604+4304 


124 


< 1.0 


0.57 


± 


0.11 


0.40 


± 


0.14 


0.17 


± 


0.14 


0.42 


± 


0.11 


CL 1604+4321 


150 


< 1.0 


0.39 


± 


0.10 


0.27 


± 


0.14 


0.12 


± 


0.14 


0.60 


± 


0.10 


RDCS 0910+5422 


146 


< 1.0 


0.38 


± 


0.10 


0.32 


± 


0.14 


0.06 


± 


0.14 


0.61 


± 


0.10 


RDCS 1252-2927 


67 


< 1.0 


0.68 


± 


0.14 


0.47 


± 


0.16 


0.21 


± 


0.16 


0.32 


± 


0.14 


RXJ 0848+4452 


214 


< 1.0 


0.30 


± 


0.09 


0.29 


± 


0.13 


0.01 


± 


0.13 


0.70 


± 


0.09 



Errors are the quadrature sum of counting statistics and classi- 
fication uncertainty. See text for details. 



APPENDIX 

COMPARING PROJECTED DENSITY ESTIMATION TECHNIQUES 

We demonstrate here that the Nearest Neighbor (NN) and Friends-of-Friends (FoF) algorithms for estimating local 
projected density yield consistent results. We also show that it is possible to construct a composite MDR or MRR 
using a combination of galaxy samples selected using spectroscopic redshifts, photometric redshifts, and samples with 
only statistically subtracted background corrections providing each of the samples are confined to specific density 
ranges. 

The FoF algorithm provides an alternative to the NN algorithm as a means to measure the morphological composition 
as a function of density. The FoF algorithm, first developed by Huchra & Geller (1982), has been used widely thereafter 
for automated detection of overdensities in galaxy catalogs. This approach has the advantage that it is not tied to a 
specific choice of N nearest neighbors although it is not the method typically used by other workers in the analysis of 
the MDR. The FoF algorithm locates groups of galaxies by identifying all the neighboring systems that lie within a 
given "percolation" length of a given galaxy. Each neighbor is then searched for galaxies within the same percolation 
length. This linking process is continued until no additional members can be identified. The percolation length and the 
local projected density are related as Iperc oc For our comparison, we use percolation lengths that correspond 

to overdensities in the range 15 < S < 1500 galaxies Mpc~^. At each overdensity, the morphological populations 
within all identified groups are tallied, weighted by the inverse of the selection function. The difference between the 
cumulative E, SO, and Sp+Irr population at two overdensities E and S — then defines the morphological populations 
at projected density S. As with the nearest neighbor algorithm, we correct the FoF densities to correspond to our 
fiducial l umin osity limit using equation |2| 

Figure lXTBl shows a comparison between the NN and FoF density estimators for the clusters for which we have reliable 
photometric redshifts (MS1054-0321, RXJ0152-1357, and RDCS1252-2927). The errorbars in all these figures include 
the uncertainties from counting statistics and from the intrinsic classification uncertainties quantified in section |31(e. (7., 
±0.06 in /e+so and /sp+irr; ±0.11 in fso and /e). The two density estimators produce consistent results over the full 
range of densities analyzed. 

Using Composite Cluster Samples 

Because not all clusters in our sample have sufficient data to compute the local density based on either spectroscopic 
or on photometric redshift values, it is also important to establish that our results are not very sensitive to the details 
used to derive the MDR {e.g., whether the sample is selected based on photometric redshift value, spectroscopic 
redshift value, or a flux- limited sample with statistical background subtraction). If we can demonstrate this then we 
can construct a composite MDR by combining for each individual cluster the sample that yields the most reliable 
available estimate {e.g., the samples with the least contamination from foreground/background objects). The clusters 
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Fig. A16.— The composite morphology - density relation for the clusters MS1054-0321, RXJ0152-1357, and RDCS1252-2927. These 
results are based on the photo-z selected samples that include galaxies with Zpf^ values in the range {z^i — Zpfi\/(1 + Zcl) ^ '^'^ph- Results 
from the A'^th nearest-neighbor (NN) and friends-of-friends (FoF) density estimators are both shown here. The density values have been 
offset by a small amount from one another for clarity. The total number of galaxies in each density bin is shown in the top panel of the 
plot. The low-density population fractions from the SDSS (Goto et al. 2003a) are denoted by the triangles along the y-axis. 



RXJ0152-1357 and MS1054-0321 allow us to perform this test as they have sufRcient information to independently 
generate the MDR from spectroscopic data, from photometric redshift data, and from flux-limited samples with a 
statistically-subtracted background correction applied. The results are shown in Figure [^TtI For densities greater 
than 80 galaxies Mpc~^ there are no significant systematic differences in any of the derived population fractions 
as a function of density. Below 80 galaxies Mpc~^, the statistically subtracted background corrected density and 
population fraction estimates become less reliable as fluctuations in the background galaxy surface density become 
comparable with the projected density (see discussion in section^. The spectroscopic samples exhibit a lower elliptical 
galaxy fraction at densities below about 50 galaxies Mpc~^ but the differences lie within the la uncertainties. The 
spectroscopic data do not provide good samphng of the very highest densities (S > 1000 galaxies Mpc~^) because 
of the number of slit masks used. However, in this regime, the photo-z and/or statistically subtracted background 
results are very reliable. We conclude that the MDR and MRR derived from a composite sample will be reliable if 
spectroscopically selected samples are limited to regions with S < 1000 galaxies Mpc^^ and samples with statistically 
subtracted background corrections are limited to regions with S > 80 galaxies Mpc~^. 
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Fig. A17. — The dependence of MDR on sample selection is demonstrated for MS1054-0321 and RXJ0152-1357. We show the population 
fractions as a function of projected density based on the samples of confirmed spectroscopic members (Spec), the samples of photo-z selected 
members (BPZ), and on the full flux-limited samples with statistically subtracted background corrections applied (SSub). The latter results 
are only shown when E > 40 galaxies Mpc~^. Errorbars include the uncertainties in counting statistics and morphological classification. 



COMPUTING MORPHOLOGICAL POPULATION FRACTIONS 

Our morphological population fractions are corrected for both contamination and incompleteness in a manner that 
depends on the sample being analyzed. For samples using only confirmed spectroscopic redshifts, the only correction 
applied is the weighting by the inverse of the redshift selection function (cf. equation^). The redshift selection function 
is computed empirically by measuring the ratio of the number of redshifts acquired to the total number of galaxies in 
magnitude bins 0.5 mag wide. The redshift selection function is dependent on the galaxy magnitude and to a lesser 
degree on the galaxy color. The latter dependence translates to a dependence on morphology for those galaxies in 
the cluster. By measuring the redshift completeness as a function of color and magnitude, we can make reasonable 
corrections for the observational selection effects. The clusters for which we have a sufficient number of redshifts to 
compute the MDR solely from a spectroscopic sample are MS1054-0321 (with 143 confirmed members, all within the 
boundaries of our ACS mosaic) and RXJ0152-1357 (with 102 confirmed members, 93 of which lie within the boundaries 
of our ACS mosaic). The redshif t sele ction functions for these two clusters is shown as a function of 1775 magnitude 
and morphological type in Figure IB 181 

For photometric redshift selected samples (MS1054-0321, RXJ0152-1357, and RDCS1252-2927), we correct the 
population fractions as follows: 

j^corr ^ J^^cor _^ j^jnissed _ j^contani ^g^^ 
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Fig. B18. — The redshift selection functions for E+SO and Sp+Irr galaxies vs. 1775 magnitude for MS1054-0321 and RXJ0152-1357. The 
larger difference between E+SO and Sp+Irr selection efficiency in RXJ0152-1357 is caused by the color selection criterion used in target 
selection for that cluster. 



where N^"-^""^ is the observed number of galaxies of morphological type T, ^v™*^'*'^'' is an estimate of the number of 
cluster members of type T that have been excluded by our photo-z selection criteria, and _^r™"*a™ ig an estimate of 
the number of non-cluster members of type T that have been included in the observed count. The number of missed 
members is computed as follows: 

jymissed _ jy^^^^^ y ftot,miss X fmiss,T (B2) 

where Ntotai is the total number of galaxies counted at a given density, ftot,miss is the fraction of cluster members 
excluded by our photo-z selection limits, and fmiss,T is the fraction of those excluded galaxies that have morphological 
type T. The value of ftot,miss is estimated by counting the number of spectroscopically confirmed members that 
have photo-z's outside our photo-z selection limits. For MS1054-0321 and RXJ0152-1357, ftot,miss = 0.38 and 0.20, 
respectively. For RDCS1252-2927, ftot.miss = 0.04. The morphological dependence of the number of true cluster 
galaxies missed by our ph oto-z selection is weak for MS1054-0321 and RDCS1252-2927 but is more significant for 
RXJ0152-1357 (see Figure IbT8)| . For MS1054-0321 and RDCS1252-2927, we set fmiss,T equal to the initial uncor- 
rected population fraction for galaxies of type T at the given density. For RXJ0152-1357 we assign fmiss,T using a 
morphological distribution that is slightly more heavily weighted to late-type galaxies at a given density to compensate 
for the spectroscopic selection bias in this cluster. Fortunately, the results are not strongly dependent on the precise 
prescriptions for setting the fmiss.T values. The number of non-cluster members contaminating our photo-z selected 
counts is computed as follows: 

^contam ^ ^^^^^^ ^ f^^^ 
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where ftot,contam is the fraction of non-cluster members included by our photo-z selection limits and fcontam,T is the 
fraction of these that have morphological type T. The value of ftot,contam is estimated by counting the number of 
spectroscopically confirmed non-cluster members that have photo-z's lying within our photo-z selection limits. For 
MS1054-0321 and RXJ0152-1357, ftot,contam = 0.09 and 0.12, respectively. For RDCS1252-2927, ftot,contam = 0.33. 
These non-cluster galaxies have a morphological distribution that is representative of the field galaxy population and 
we thus set fcontam,E = 0.10, fcontam.so = 0.25, and fcontam,Sp+irr = 0.65. Notc that bccausc we only have 35 
spectroscopic redshifts for RDCS1252-2927, the incompleteness and contamination estimates have larger uncertainties 
than those in MS1054-0321 or RXJ0152-1357. 

For population fractions derived from magnitude limited samples using only statistically-subtracted background 
counts (e.^., as in RDCS0910-f 5422), the morphological distribution of the background is assumed to be representative 
of the field galaxy population and we thus divide the expected background counts into E, SO, and Sp-|-Irr components 
assuming the same fractions given above: — 0.10, fso — 0.25, &ndfsp+irr — 0.65. 
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